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Abstract  Studying  the  evolution  of  magnetic  clouds  entrained  in  coronal  mass 
ejections  using  in-situ  data  is  a  difficult  task  since  only  a  limited  number  of 
observational  points  is  available  at  large  heliocentric  distances.  Remote  sensing 
observations  can,  however,  provide  important  information  for  events  close  to 
the  Sun.  In  this  work  we  estimate  the  flux  rope  orientation  first  in  the  close 
vicinity  of  the  Sun  (2-20  Rq)  using  forward  modeling  of  STEREO/SECCHI  and 
SOHO/LASCO  coronagraph  images  of  coronal  mass  ejections  and  then  in-situ 
using  Grad-Shafranov  reconstruction  of  the  magnetic  cloud.  Thus,  we  are  able  to 
measure  changes  in  the  orientation  of  the  erupted  flux  ropes  as  they  propagate 
from  the  Sun  to  1  AU.  We  present  both  techniques  and  use  them  to  study 
15  magnetic  clouds  observed  during  the  minimum  following  Solar  Cycle  23  and 
the  rise  of  Solar  Cycle  24.  This  is  the  first  multievent  study  to  compare  the  three- 
dimensional  parameters  of  CMEs  from  imaging  and  in-situ  reconstructions.  The 
results  of  our  analysis  confirm  earlier  studies  showing  that  the  flux  ropes  tend  to 
deflect  towards  the  solar  equatorial  plane.  We  also  find  evidence  of  rotation  on 
their  travel  from  the  Sun  to  1  AU.  In  contrast  to  past  studies,  our  method  allows 
one  to  deduce  the  evolution  of  the  three-dimensional  orientation  of  individual 
flux  ropes  rather  than  on  a  statistical  basis. 
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etary;  Magnetic  fields,  Models 


1.  Introduction 

Coronal  mass  ejections  (CMEs)  are  massive  bursts  of  plasma  and  magnetic  field 
from  the  Sun  into  the  interplanetary  space.  Interplanetary  coronal  mass  ejections 
(ICMEs)  are  the  heliospheric  counterparts  of  the  CMEs.  ICMEs  are  one  the 
main  drivers  of  space  weather  (Tsurutani  et  al. ,  1988,  Huttunen  et  al,  2002, 
Zhang  et  al.,  2007).  ICMEs  show  a  variety  of  signatures  in  in-situ  observations 
at  1  AU.  Some  of  them  such  as  the  enhancement  of  the  magnetic  field,  smooth 
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monotonic  rotation  of  the  magnetic  field  through  a  large  angle,  low  proton 
temperature  and  low  plasma  /3,  indicate  the  existence  of  a  flux  rope  structure 
within  the  body  of  the  ICME.  ICMEs  with  embedded  magnetic  flux  ropes  are 
called  magnetic  clouds  (MCs).  The  list  of  signatures  found  in  in- situ  measure¬ 
ments  of  MCs  can  be  found  in  the  paper  of  Zurbuchen  and  Richardson  (2006). 
Approximately  one-third  of  all  ICMEs  observed  at  1  AU  show  magnetic  cloud 
signatures  (Gosling,  1990),  but  there  are  indirect  evidences  that  all  ICMEs  might 
have  central  flux  ropes.  Krall  (2007)  conluded  after  analyzing  Solar  Maximum 
Mission  (SMM)  coronagraph  data  that  all  CMEs  in  his  sample  may  share  a 
similar  hollow-flux-rope  structure  (see  also  Vourlidas  et  al. ,  2013,  this  volume). 
As  suggested  by  Jian  et  al.  (2006),  the  absence  of  MC  signatures  in  many  ICMEs 
may  be  a  positional  effect:  their  analysis  showed  that  in  about  two-thirds  of  the 
cases  the  spacecraft  encounters  the  ICME  so  far  from  the  center  that  the  central 
flux  rope  is  not  identifiable.  Since  the  southward  magnetic  field  has  the  strongest 
influence  on  the  Earth’s  magnetosphere,  the  geoeffectivity  of  MCs  depends  on 
their  orientation.  The  orientation  of  the  MC  in  the  low  solar  corona  and  at  1  AU 
can  differ  greatly,  so  it  is  crucial  to  understand  how  the  orientation  of  MCs 
changes  on  their  journey  from  the  Sun  to  1  AU.  This  knowledge  will  help  us  to 
better  understand  the  processes  affecting  CME  propagation  in  the  heliosphere, 
improve  the  solar  and  heliosphere  magnetohydrodynamics  simulation  techniques 
and  hence  improve  space  weather  forecasting. 

The  orientation  of  MCs  during  their  journey  from  the  Sun  to  1  AU  can  change 
significantly.  The  change  of  the  orientation  can  be  decomposed  into  latitudi¬ 
nal  deflection,  longitudinal  deflection  and  rotation.  The  latitudinal  deflection 
of  CMEs  in  the  low  corona  was  reported  already  by  MacQueen  et  al.  (1986). 
According  to  their  work,  CMEs  tend  to  deflect  towards  the  solar  equatorial  plane 
by  about  2  degrees  on  average  while  travelling  from  2 Rq  to  4i?Q  during  solar 
minimum.  Plunkett  et  al.  (2001)  found  that  the  latitude  distribution  of  CMEs 
in  the  outer  corona  near  solar  minimum  is  very  different  to  the  distribution  in 
the  inner  corona,  suggesting  that  the  propagation  of  CMEs  in  the  inner  corona 
is  controlled  by  the  large-scale  solar  magnetic  field  which  tends  to  push  the 
CME  towards  the  equatorial  streamer  belt  as  it  propagates  outward.  On  the 
other  hand,  CMEs  show  no  average  latitudinal  deflection  during  solar  maximum. 
Cremades  et  al.  (2006)  has  shown  that  the  deviation  of  CMEs  with  respect  to 
their  source  regions  is  always  equatorward  near  solar  minimum,  while  deviations 
to  higher  latitudes  are  also  frequent  during  solar  maximum.  They  also  found  a 
significant  correlation  of  the  deviation  with  the  number  of  coronal  holes,  their 
area  and  their  distance  to  the  CME  source  regions.  Kilpua  et  al.  (2009)  pointed 
that  the  spatial  dimensions  of  the  erupting  CME  might  play  a  role  in  determining 
whether  a  CME  will  be  deflected  towards  the  equator,  i.e.  slower  CMEs  with 
wide  longitudinal  extent  could  not  penetrate  through  the  background  coronal 
fields. 

The  evidence  for  longitudinal  deflection  is  based  on  the  statistics  of  the  solar 
source  regions  of  geoeffective  halo  CMEs  (i.e.  CMEs  that  arrived  to  the  Earth 
and  produced  geomagnetic  storms).  It  was  found  by  Wang  et  al.  (2002)  that 
the  East- West  distribution  of  the  CME  sources  is  asymmetrical  -  the  number 
of  geoeffective  halo  CMEs  originating  from  the  western  hemisphere  is  larger 
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by  57%  compared  to  the  ones  originating  from  the  eastern  hemisphere.  In  the 
paper  by  Wang  et  al.  (2011)  the  CME  deflections  were  classified  into  three 
types:  (i)  asymmetrical  expansion  of  the  CME,  (ii)  nonradial  ejection  and  (iii) 
deflected  propagation.  Deflected  propagation  in  the  lower  corona  can  be  caused 
by  the  interaction  of  the  CME  with  other  neighbouring  large-scale  magnetic  field 
structures  such  as  coronal  holes. 

MCs  can  also  experience  rotation  on  their  journey  to  1  AU.  The  possible 
reason  for  MC  rotation  can  be  the  interaction  with  neighbouring  magnetic 
field  structures  and  the  kink  instability  (Torok  and  Kliem,  2003).  But  CME 
rotation  is  not  expected  at  large  distances  when  the  ambient  magnetic  field 
is  weak  (Lynch  et  al.,  2009).  Sometimes  MCs  can  suffer  rather  rapid  rotation, 
however.  Vourlidas  et  al.  (2011)  reported  a  MC  rotating  at  a  rate  of  60°  per  day 
in  the  low  corona.  Based  on  statistical  evidence  MCs  seem  to  rotate  towards  the 
heliospheric  current  sheet  (HCS)  so  that  the  MC  stays  aligned  with  the  local 
HCS  as  shown  by  Yurchyshyn  (2008)  and  Yurchyshyn  et  al.  (2009). 

It  is  difficult  to  study  the  evolution  of  MCs  from  the  Sun  to  1  AU,  since  there 
is  only  a  limited  number  of  observing  points  for  analysis  of  such  a  complicated 
and  large  structure.  The  launch  of  the  STEREO  mission  (Kaiser  et  al.,  2008) 
made  it  possible  to  obtain  stereoscopic  coronagraph  images  of  the  Sun  which 
gave  rise  to  forward  modeling  (Thcrnisien  et  al.,  2009,  2011).  Paired  with  in-situ 
measurements  it  allows  to  obtain  the  geometrical  parameters  of  MCs  in  different 
temporal  and  spatial  stages  of  their  evolution.  In  this  work  we  study  the  evolution 
of  MCs  registered  during  the  minimum  following  Solar  Cycle  23  and  the  rise  of 
Solar  Cycle  24.  We  use  forward  modeling  technique  to  study  the  properties  of 
MCs  close  to  the  Sun  (5-20  Rq)  and  Grad-Shafranov  reconstruction  at  1  AU. 

Our  paper  is  organized  as  follows:  in  Section  2  we  will  describe  our  method  of 
estimating  the  orientation  of  MCs  using  one  example  from  our  list  of  15  events, 
in  Section  3  we  present  the  results  of  our  analysis  and  discuss  them  in  Section 
4. 

2.  Methodology  using  an  example  event 

A  coronagraph  is  a  telescope  pointed  at  the  Sun  with  an  occulter  blocking  the 
solar  disk  light.  The  solar  disk  is  about  106  —  109  times  brighter  than  the  inner 
corona,  thus  the  occulter  is  necessary  to  observe  the  solar  corona  with  sufficient 
contrast  to  reveal  faint  structures.  There  are  three  coronagraphs  on  board  the 
SOHO  spacecraft  with  fields  of  view  of  1.1-3  Rq  (Cl),  2.2-6  Rq  (C2)  and  3.5- 
30  Rq  (C3),  part  of  the  LASCO  ( Large  Angle  and  Spectrometric  Coronagraph) 
experiment  (Brueckner  et  al.,  1995).  The  STEREO  spacecraft  are  supplied  with 
two  coronagraphs  each  with  fields  of  view  of  1.5-4  Rq  (COR1)  and  2.5-15  Rq 
(COR2)  and  a  heliospheric  imager  (HI)  pointed  towards  the  Sun-Earth  line,  part 
of  the  SECCHI  ( Sun  Earth  Connection  Coronal  and  Heliospheric  Investigation : 
Howard  et  al.  (2008))  package.  In  our  study  we  tracked  each  CME  from  2  to 
20-40  Rq  using  LASCO  C2  and  C3  and  SECCHI  COR2  and  HI  telescopes. 
We  required  that  each  CME  was  captured  in-situ  by  at  least  one  spacecraft 
(STEREO  and/or  Wind). 


SOLA:  Isavnin_etal_2012a.tex;  5  December  2012;  4:31;  p.  3 


A.  Isavnin,  et  al. 


Figure  1.  Scheme  of  the  event  suitable  for  our  analysis. 


The  method  we  developed  in  this  paper  to  study  the  evolution  of  a  MC  is  based 
on  the  forward  modeling  of  coronagraph  data  (Thernisien,  Vourlidas,  and  Howard,  2009)| 
and  Grad-Shafranov  reconstruction  of  in-situ  MCs  (Hu  and  Sonnerup,  2002). 

We  present  this  method  alongside  an  example  event  of  the  CME  on  4  Novem¬ 
ber  2010  and  the  associated  in-situ  MC  measured  on  8  November  2010  by  the 
STEREO-B  spacecraft.  The  schematic  representation  of  this  event  is  depicted 
on  Figure  1. 

We  assume  that  the  invariant  axis  of  a  MC  lies  in  a  single  plane  (Figure  2), 
so  that  when  the  MC  changes  its  orientation,  the  plane  containing  the  invariant 
axis  changes  its  orientation  accordingly.  We  will  refer  to  it  as  the  MC  plane  from 
now  on.  The  orientation  of  the  MC  can  be  characterised  by  the  normal  to  the 
MC  plane  and  the  direction  from  the  Sun  to  the  apex  of  the  MC  (Figure  2). 

The  multi-spacecraft  forward  modeling  (FM)  was  introduced  by  Thernisien,  Vourlidas,  and  Howard| 
(2009)  and  is  implemented  in  the  SolarSoft  package.  The  technique  is  based  on 
the  fitting  of  a  three-dimensional  hollow-croissant-shaped  structure  to  stereo¬ 
scopic  coronagraph  images  of  a  CME.  It  is  preferable  to  use  all  three  spacecraft 
observations  for  a  fitting  procedure  since  at  least  one  of  them  will  observe  the 
CME  as  halo  or  partial  halo  resulting  in  faint  coronagraph  images.  Figure  3 
shows  the  FM  fitting  for  the  CME  on  4  November  2010.  The  separation  between 
the  STEREO  spacecraft  was  about  166°,  so  both  STEREO  spacecraft  observed 
partial  halo  CME  (backside  partial  halo  CME  for  STEREO-A).  LASCO,  at  the 
same  time,  observed  the  CME  on  the  Eastern  part  of  the  limb.  The  output 
of  FM  that  we  are  interested  in  are  the  direction  of  radial  CME  propagation 
(0fm,  <^fm)  and  its  rotation  angle  7fm-  These  outputs  are  given  in  the  Stonyhurst 
coordinate  system  (Thompson,  2006).  The  direction  of  CME  propagation  ?fm 
and  the  normal  to  the  MC  plane  £lfm  in  the  lower  corona  are  determined  then 
as 


f  FM  =  COS0FM  COS^FMex  +  COS0FM  sin^FMGy  +  sin0FMGz,  (1) 

nFM  =  rotated ?fm  x  ey]  x  rFM, ?fm, Tfm),  (2) 

where  we  defined  vector  rotation  operator  rotate  as 

rotate (v,  a,  7)  =  v  cos  7  +  (v  ■  a)  (1  —  cos  7)a  +  [ax  v]  sin  7,  (3) 
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Figure  2.  Three-dimensional  scheme  of  magnetic  cloud.  ?pM  is  the  global  axis  of  the  flux 
rope,  nFR  is  the  local  invariant  axis  of  the  flux  rope  at  the  point  of  intersection  with  the 
spacecraft,  Pgsr  's  the  impact  parameter. 


which  rotates  v  around  a  by  angle  7  counterclockwise. 

After  having  determined  the  orientation  of  the  CME  close  to  the  Sun  we 
turn  our  attention  to  in-situ  observations  near  1  AU.  The  Grad-Shafranov  re¬ 
construction  (GSR)  (Hu  and  Sonnerup,  2002)  is  used  for  the  estimation  of  the 
local  direction  of  the  invariant  axis  of  MC  and  reconstruction  of  a  slice  of  the 
MC.  We  use  the  modified  version  of  GSR  described  in  Isavnin  et  al. ,  (2011).  The 
output  we  are  interested  in  are  the  local  direction  of  the  invariant  axis  of  the 
flux  rope  (0GSR)  ^GSr) 

nFR  =  cos  6>gsr  cos  </>gsr  ex  +  cos  0gsr  sin  <j>  GSR^y  +  sin  $GSRez ,  (4) 

and  the  impact  parameter  pgsr-  The  impact  parameter  is  the  measure  of  the 
closest  approach  of  the  spacecraft  to  the  invariant  axis  of  the  flux  rope,  calculated 
either  as  the  distance  of  the  closest  approach  (in  astronomical  units)  or  the 
distance  of  the  closest  approach  divided  by  the  radius  of  the  flux  rope  cross- 
section.  Knowing  the  position  of  the  spacecraft  which  crossed  the  MC  (0sc>  </>sc) 

r  sc  =  cos  0sc  cos  0scex  +  cos  9Sc  sin  </>scey  +  sin  0scez  (5) 

and  the  impact  parameter  pgsr  we  can  estimate  the  vector  pointing  to  the 
part  of  the  of  flux  rope  closest  to  the  spacecraft  trajectory  using  the  following 
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Figure  3.  Example  of  forward  modeling  (lower  panels)  of  the  CME  of  4  November  2010 
(upper  panels).  The  gray-scale  images  from  left  to  right  are  the  coronagraph  images  obtained  by 
STEREO-B/SECCHI  COR2,  SOHO/LASCO  C3  and  STEREO-A/SECCHI  COR2  telescopes. 


equation: 

.  .  ~  o  •  sign(hFR- ey)pGSR\  ,a, 

rFR  =  rotate  I  rSc,  nFR,  2  arcsm - -  I  .  (6) 

The  vector  defined  by  Equation  (6)  lies  in  the  MC  plane.  We  can  estimate  the 
normal  to  the  MC  plane  at  1  AU  as 


nGSR  =  nFR  x  fFR.  (7) 

The  MC  associated  with  the  CME  in  example  event  was  registered  in  situ 
by  the  STEREO-B  spacecraft  on  8  November  2010.  Figure  4  shows  the  residual 
map  and  the  reconstructed  magnetic  field  map  for  this  event.  The  residual  map 
shows  the  process  of  the  search  for  the  local  direction  of  the  invariant  axis  of 
the  flux  rope.  It  represents  the  hemisphere  of  all  possible  orientations  of  the  axis 
and  the  direction  with  the  minimal  residue  is  the  estimated  invariant  axis  of  the 
flux  rope.  The  reconstructed  magnetic  field  map  is  essentially  the  cross-section 
of  the  flux  in  the  vicinity  of  the  spacecraft  trajectory.  Black  arrows  show  the 
magnetic  field  measured  in  situ  and  projected  onto  the  plane  perpendicular  to 
the  invariant  axis  of  the  flux  rope.  Black  contour  lines  denote  the  equipotential 
levels,  where  the  absolute  values  of  the  vector  potential  are  considered.  The 
white  dot  represents  the  invariant  axis  of  the  flux  rope.  The  thick  white  contour 
line  shows  the  boundary  of  the  unperturbed  part  of  the  flux  rope. 
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Figure  4.  Residual  map  (left)  and  magnetic  field  map  (right)  for  the  MC  registered  on  8 
November  2010  by  STEREO-B  spacecraft.  The  projected  RTN  coordinate  axes  are  Rrtn 
(cyan),  Trtn  (magenta)  and  Nrtn  (yellow) 


low  pressure  high  pressure 

Figure  5.  Interaction  of  a  slow  (left)  and  fast  (right)  CME  with  Parker-spiral-structured  solar 
wind  (adapted  from  Wang  et  al .,  2004). 


Wang  et  al.  (2004)  introduced  the  kinetic  interpretation  of  the  longitudinal 
deflection  of  MCs  which  conforms  with  the  statistics  of  observations  of  geoeffec- 
tive  MCs  and  their  sources  (Wang  et  al.,  2002).  According  to  these  authors  the 
longitudinal  deflection  of  MCs  is  caused  by  the  interaction  of  the  MC  with  the 
Parker-spiral-structured  solar  wind  (see  Figure  5).  If  an  MC  propagates  slower 
than  the  background  solar  wind  it  is  pushed  by  the  faster  solar  wind  thus  getting 
the  westward  component  of  the  force.  On  the  other  hand,  a  fast  MC  is  blocked 
by  the  slower  background  solar  wind  thus  getting  the  eastward  component  of 
the  force.  Longitudinal  deflection  can  be  estimated  using  the  following  equation: 

A^nfe-vfc)'1AU-  (8> 

where  Vmc  is  the  average  speed  of  the  MC  propagation,  Vsw  is  the  average 
velocity  of  the  background  solar  wind  and  f 1  «  2.7  x  10_6rads_1  is  the  angular 
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velocity  of  the  Sun’s  rotation.  The  positive  value  of  A <j>  represents  westward 
deflection,  while  the  negative  value  of  A cj>  shows  eastward  deflection  of  the  MC. 
MCs  which  propagate  with  the  velocity  close  to  the  background  solar  wind 
velocity  do  not  experience  longitudinal  deflection  on  average.  It  should  be  noted 
that  this  is  a  rough  estimate  since  the  speed  of  the  MC  can  change  on  its  journey 
from  the  Sun  to  1  AU.  Since  we  are  studying  only  events  registered  near  the  solar 
minimum,  most  of  which  are  slow  CMEs,  the  simplicity  of  Equation  (8)  should 
not  affect  the  results  of  the  analysis  significantly.  For  our  example  event  the  MC 
propagation  velocity  at  1  AU  is  slightly  lower  than  the  background  solar  wind. 
So  we  estimate  the  westward  deflection  using  Equation  (8)  to  be  Acj)  =  1°. 

Longitudinal  deflection  is  taken  into  account  by  rotating  the  initial  MC  di¬ 
rection  fpM  by  Acj)  around  e~: 

?fm  =  rotate( ?Fr,  ez,  Acj>).  (9) 

Since  fFM  is  directed  along  the  global  axis  of  the  MC  through  its  apex,  the 
possible  rotation  of  the  MC  around  its  axis  will  not  affect  fpM.  Then  only 
latitudinal  deflection  can  change  r'FM: 

f fm  =  rotate(r'FM,e'y,  Ad),  (10) 

where  =  ez  x  rFM  and  A 0  is  latitudinal  deflection  angle.  The  vectors  ?pM 
and  fpR  both  he  in  the  MC  plane  (see  Figure  2)  and  its  normal  at  1  AU  can  be 
estimated  as 


**FM  —  fFR  X  ?FM-  (11) 

Latitudinal  deflection  is  determined  as  such  an  angle  Ad  for  which  ncsR  •  fipM  = 
1.  The  positive  value  of  A 9  represents  North-to-South  deflection  and  negative 
value  shows  South-to-North  rotation.  For  our  example  event  we  found  the  lati¬ 
tudinal  deflection  to  be  A 9  =  —17°.  The  MC  deflected  from  the  Southern  to  the 
Northern  hemisphere,  crossing  the  helioequatorial  plane. 

After  estimating  Acj)  and  A 9  deflections  the  rotation  can  be  calculated  using 
the  following  equation: 


Ay  =  arctan  (  ^FR  f z  )  -  7fm,  (12) 

VrFR  •  e’yJ 

where  e'z  =  fpM  x  e^.  The  positive  value  of  Ay  represents  counterclockwise 
rotation,  while  the  negative  value  of  Ay  shows  clockwise  rotation  of  the  MC 
around  its  global  axis.  Using  Equation  (12)  we  estimate  a  rotation  by  Ay  =  30° 
for  our  example  event.  It  should  be  noted  that  the  FM  fit  contains  no  information 
about  polarity  or  chirality  of  the  flux  rope.  Thus,  the  estimated  initial  rotation 
angle  of  the  flux  rope  has  an  ambiguity  of  180  degrees.  This  implies  that  Equation 
(12)  shows  the  change  of  the  tilt  angle  of  the  MC  plane,  but  not  the  real  direction 
or  full  amount  of  rotation. 
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Figure  6.  Diagrams  of  deflection  (left)  and  rotation  (right)  of  the  studied  MCs.  Red  squares 
denote  the  direction  of  MC  propagation  (left)  and  the  direction  of  the  normal  to  the  MC  plane 
(right)  close  to  the  Sun,  i.e.  in  the  initial  position.  Green  diamonds  denote  the  direction  of 
MC  propagation  (left)  and  the  direction  of  the  normal  to  the  MC  plane  (right)  at  1  AU,  i.e. 
in  the  final  position. 


3.  Results 

For  the  analysis  we  selected  the  MCs  which  were  observed  during  the  minimum 
following  the  Solar  Cycle  23  and  the  rise  of  the  Solar  Cycle  24,  i.e.  years 
2008-2010.  Such  a  choice  is  based  on  several  considerations,  i.e.  sufficiently 
large  angular  distance  (at  least  30-40  degrees)  between  the  STEREO  spacecraft 
which  was  launched  in  2006  is  required  for  the  FM  technique  results  to  be 
reliable  (Thernisien,  Vourlidas,  and  Howard,  2009).  The  separation  between  the 
STEREO  and  SOHO  spacecraft  grew  from  26°  to  87°  during  that  period.  The 
minimum  of  solar  activity  is  also  characterised  by  a  less  dynamic  structure  of 
the  heliospheric  current  sheet  (HCS)  and  coronal  holes  which  facilitate  studies 
of  interaction  of  an  MC  with  these  structures. 

Similar  to  other  flux  rope  fitting  techniques,  GSR  works  best  for  small  impact 
parameter  events  (Isavnin  et  ai,  2011),  so  we  conducted  analysis  only  for  MCs 
which  were  crossed  close  to  their  invariant  axis  by  the  spacecraft  at  1  AU.  We 
have  selected  15  events  for  our  analysis,  the  results  of  analysis  are  summarized 
in  Table  1. 

The  diagrams  in  Figure  6  visualize  the  estimated  deflection  and  rotation 
experienced  by  the  analyzed  MCs  when  they  propagated  from  the  Sun  to  1  AU. 
Also,  the  last  three  columns  of  Table  1  give  the  longitudinal  and  latitudinal 
deflection  angles  (calculated  from  equations  (8)-(ll))  and  the  rotation  angle 
estimated  from  Equation  (12)  for  each  event.  In  the  left  diagram  of  Figure  6 
the  ^-components  of  the  lines  connecting  the  red  squares  and  green  diamonds 
show  the  amount  of  longitudinal  deflection,  while  the  y-components  of  the  lines 
represent  the  amount  of  latitudinal  deflection.  In  the  right  panel,  the  longer  the 
angular  distance  covered  by  a  curve  between  the  red  squares  and  green  diamonds 
the  more  the  MC  has  rotated  during  its  interplanetary  propagation. 
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As  shown  by  Figure  6  and  Table  1  for  all  studied  events  the  longitudinal 
deflection  was  small,  less  than  6  degrees.  The  latitudinal  deflection  was  clearly 
larger  than  the  longitudinal,  exceeding  10°  for  the  majority  of  events  with 
the  largest  deflection  of  35°  (event  #10).  According  to  our  study  14  out  of 
15  analyzed  MCs  deflected  towards  the  helioequatorial  plane,  often  crossing 
it.  The  only  exception  is  event  #6  for  which  the  deflection  by  2°  away  from 
helioequatorial  plane  was  found.  This  event  is  also  the  one  with  the  smallest 
amount  of  latitudinal  deflection  among  the  15  events  we  studied. 

The  amount  of  rotation  of  the  studied  MCs  ranged  from  events  that  experi¬ 
enced  practically  no  rotation  (event  #12)  to  events  that  rotated  as  much  as  80 
degrees  (event  #4).  The  average  absolute  value  of  rotation  angle  for  our  data 
set  was  31  degrees. 

From  the  Figure  6(left)  it  can  be  seen  that  the  directions  of  propagation  of 
14  out  of  15  MCs  deflected  to  or  stayed  in  the  latitude  range  of  dstonyhurst  = 
[— 12°,12°]  after  their  travel  from  the  Sun  to  1  AU.  Thus,  MCs  seem  to  align 
the  direction  of  their  expansion  with  the  solar  equatorial  plane. 


4.  Discussion 

We  have  presented  the  first  detailed  study  of  the  three-dimensional  evolution 
of  the  orientation  of  MCs  from  the  Sun  to  1  AU.  The  analysis  makes  use  of 
coronagraph  images  of  CMEs  and  in-situ  observations  of  erupted  flux  ropes  and 
is  capable  of  more  precise  identification  of  the  three-dimensional  orientation  of 
individual  flux  ropes  and  their  evolution  than  was  possible  in  the  past.  This 
method  utilizes  the  FM  and  GSR  techniques  for  estimation  of  the  flux  rope 
orientation  close  to  the  Sun  and  at  1  AU,  respectively.  The  approach  is  con¬ 
strained  by  certain  limitations  on  the  quality  of  the  studied  events.  The  flux 
rope  structure  within  the  CME  associated  with  the  analyzed  MC  has  to  be  clear 
enough  to  be  easily  distinguishable  in  the  coronagraph  images.  The  separation 
between  the  STEREO  spacecraft  has  to  be  large  enough  (at  least  30-40  degrees) 
to  get  three-dimensional  representation  of  the  CME  and  hence  a  better  fit  of  FM. 
The  impact  parameter  of  the  MC  observed  at  1  AU  has  to  be  small  to  estimate 
the  orientation  of  MC  at  1  AU  more  precisely  using  GSR  technique,  which  works 
better  for  small  impact  parameters. 

Based  on  these  considerations  we  have  selected  15  clear  events  during  the  years 
2008-2010  and  used  our  method  to  study  their  three-dimensional  evolution  from 
the  Sun  to  1  AU.  The  events  were  observed  during  the  minimum  following  Solar 
Cycle  23  and  the  rise  of  Solar  Cycle  24.  Our  analysis  shows  that  MCs  tend  to 
deflect  towards  the  solar  equatorial  plane  on  their  journey  from  the  Sun  to  1  AU. 
This  result  is  in  agreement  with  previous  statistical  studies  by  MacQueen  et  al. 
(1986),  Plunkett  et  al.  (2001)  and  Cremades  et  al.  (2006),  though  in  our  work 
we  were  able  to  calculate  the  evolution  of  MCs’  orientation  directly  from  the 
multi-spacecraft  observations.  In  this  study  we  have  not  considered  the  source 
regions  of  CMEs  but  we  started  tracking  MCs  from  2  Rq  ,  hence  we  showed  that 
the  latitudinal  deflection  of  MC  can  happen  not  only  in  the  lower  corona,  but 
the  orientation  of  MC  continues  to  evolve  all  the  way  to  1  AU.  A  possible  reason 
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for  the  latitudinal  deflection  may  be  the  kinematic  interaction  between  CMEs 
and  fast  solar  wind. 

The  studied  events  showed  very  little  longitudinal  deflection.  As  discussed  in 
the  Introduction,  Wang  et  al.  (2002)  observed  an  obvious  East- West  asymmetry 
in  the  source  region  distribution  of  geoeffective  halo  CMEs.  While  the  sample  of 
Wang  et  al.  (2002)  covers  the  rising  phase  of  solar  activity  and  solar  maximum, 
our  study  period  coincides  with  relatively  low  solar  activity  conditions,  and  thus 
most  of  our  CMEs  were  slow  and  embedded  into  solar  wind  of  the  speed  close 
to  the  CME  speed,  i.e.  were  not  pushed  by  fast  solar  wind  or  blocked  by  slow 
solar  wind.  Thus,  it  is  possible  that  our  data  set  represents  events  that  were  not 
much  influenced  by  the  ambient  solar  wind  flow.  We  also  expect  that  forward 
modeling  based  on  the  coronagraph  data  gives  a  more  reliable  estimate  of  the 
CME  propagation  direction  than  the  analysis  of  their  source  regions  based  on 
the  solar  disk  observation.  Our  results  are  also  consistent  with  the  analysis  by 
Rodriguez  et  al.  (2011)  who  found  that  the  predictions  of  ICME  detections  based 
on  the  forward  modeling  of  STEREO / COR2  data  matched  well  with  the  actual 
in-situ  observations.  Our  results  thus  imply  that  at  least  near  solar  minimum 
the  CME  propagation  direction  in  longitude  can  be  predicted  accurately  based 
on  the  coronagraph  data  and  this  direction  does  not  change  significantly  from  a 
few  tens  of  solar  radii  from  the  Sun  to  1  AU.  It  should  be  noted  that  although 
Equation  (8)  is  a  rough  estimate  of  the  longitudinal  deflection,  it  does  not  affect 
the  accuracy  of  the  presented  technique  dramatically.  For  instance,  introducing 
an  error  8A<f>  =  10°  to  the  longitudinal  deflection  estimate  leads  to  the  average 
error  5A9  =  3°  of  latitudinal  deflection  estimate  and  dAy  =  0.1°  of  rotation 
angle  estimate  for  the  studied  MCs. 

The  rotation  of  the  studied  MCs  could  be  caused  by  the  disconnection  of 
one  of  the  flux  rope  footpoints  early  in  the  eruption  (Vourlidas  et  al..,  2011, 
Nieves-Chinchilla  et  al..,  2012),  interaction  with  large-scale  magnetic  structures 
in  the  solar  wind  and  HCS  (Yurchyshyn,  2008) .  Analysis  of  these  possibilities  in 
relation  to  each  event  will  be  the  subject  of  our  upcoming  research. 
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Table  1.  The  results  of  the  analysis  of  MC  deflection  for  15  events.  The  CME  part  of  the  table  represents  the  results  of  FM  analysis  at  5  —  20-R©, 
MC  part  of  the  table  represents  the  results  of  GSR  analysis  at  1  AU,  the  A  orientation  represents  the  estimated  deflection  angles.  The  columns 
from  left  to  right  are  as  follows:  #  -  event  number,  date  -  date  the  CME  event  start,  <j> fm,  #fm  -  direction  of  the  CME  propagation,  7fm  -  initial 
rotation  angle  of  the  CME,  date  -  date  of  MC  registration  at  1  AU,  SC  -  spacecraft  which  observed  the  MC  at  1  AU,  </>sCj  $SC  _  coordinates 
of  the  spacecraft,  </>gsr?  #gsr  -  direction  of  the  invariant  axis  of  the  MC,  pgsr  -  impact  parameter,  A</>  -  longitudinal  deflection  (positive  for 
East- to- West),  A 6  -  latitudinal  deflection  (positive  for  North-to-South),  A 7  -  change  of  the  rotation  angle  (positive  for  counterclockwise). 


# 

date 

CME 

<t>  FM 

0FM 

7FM 

date 

SC 

</>sc 

MC 

0SC 

<^GSR 

0GSR 

PGSR  [AU] 

A  orientation 

A cf)  |  A6  |  A7 

1 

2008-06-02 

-31.85 

-2.23 

-48.49 

2008-06-06 

STB 

-24.97 

-3.77 

-152.9 

58.0 

-0.0045 

-0.58 

-13.62 

-14.20 

2 

2008-07-07 

-22.39 

-18.88 

-2.20 

2008-07-10 

STB 

-27.20 

0.10 

88.3 

-28.6 

0.0195 

-0.05 

-14.65 

-29.16 

3 

2008-08-31 

3.10 

-9.85 

-1.64 

2008-09-03 

Wind 

0.08 

7.18 

-107.4 

12.9 

0.0085 

-0.77 

-15.92 

-14.00 

4 

2008-12-12 

8.19 

7.20 

63.11 

2008-12-17 

Wind 

0.08 

-1.17 

115.7 

-15.8 

-0.0279 

-0.31 

9.26 

-80.26 

5 

2008-12-27 

-38.37 

16.37 

13.74 

2008-12-31 

STB 

-45.20 

3.40 

24.5 

-12.7 

0.0264 

1.66 

16.62 

-28.03 

6 

2009-09-27 

-52.16 

4.82 

15.49 

2009-10-02 

STB 

-56.73 

5.60 

75.1 

-42.7 

-0.0031 

-5.38 

-1.92 

-64.17 

7 

2010-01-15 

-72.40 

6.00 

15.69 

2010-01-20 

STB 

-69.07 

3.90 

-180.0 

-25.8 

-0.0045 

-0.15 

3.54 

10.54 

8 

2010-02-01 

64.54 

-17.59 

24.64 

2010-02-05 

STA 

64.10 

-6.20 

-51.9 

-31.2 

-0.0163 

-0.16 

-14.76 

10.31 

9 

2010-04-03 

4.60 

-23.66 

7.93 

2010-04-05 

Wind 

-0.24 

-6.30 

140.6 

-3.4 

0.0144 

-3.93 

-16.59 

-21.31 

10 

2010-05-23 

12.25 

6.23 

54.98 

2010-05-28 

Wind 

0.03 

-1.02 

67.6 

-60.1 

0.0007 

4.05 

34.75 

67.43 

11 

2010-05-27 

72.92 

-11.73 

62.35 

2010-05-31 

STA 

71.70 

6.50 

132.7 

0.1 

0.0165 

-2.28 

-17.46 

-65.27 

12 

2010-06-13 

97.22 

20.04 

-17.36 

2010-06-16 

STA 

73.68 

7.30 

-3.4 

15.2 

-0.0078 

-0.25 

18.85 

1.66 

13 

2010-11-04 

-76.07 

-6.56 

11.06 

2010-11-08 

STB 

-82.20 

6.40 

-139.0 

-35.3 

0.0341 

1.11 

-16.84 

29.69 

14 

2010-12-12 

48.40 

-16.87 

6.02 

2010-12-15 

STA 

85.20 

-7.30 

10.1 

10.9 

-0.0158 

2.00 

-17.86 

-21.52 

15 

2010-12-12 

-87.88 

-10.46 

-11.40 

2010-12-17 

STB 

-87.30 

7.30 

-166.0 

0.2 

0.0158 

2.78 

-2.02 

-9.85 

A.  Isavnin,  et  al. 


Three-dimensional  evolution  of  erupted  flux  ropes  from  the  Sun  (2-20  Rq)  to  1  AU 


Acknowledgements  The  work  of  A.  Isavnin  and  E.  Kilpua  was  supported  by  the  Academy 
of  Finland.  The  work  of  A.  Vourlidas  is  supported  by  NASA  contract  S-136361-Y  to  the  Naval 
Research  Laboratory.  LASCO  was  constructed  by  a  consortium  of  institutions:  NRL  (USA), 
MPI  fur  Aeronomie  (Germany),  LAS  (France)  and  University  of  Birmingham  (UK).  The 
SECCHI  data  are  produced  by  an  international  consortium  of  the  NRL,  LMSAL  and  NASA 
GSFC  (USA),  RAL  and  University  of  Birmingham  (UK),  MPS(Germany),  CSL  (Belgium), 
IOTA  and  IAS  (France). 


References 

Brueckner,  G.E.,  Howard,  R.A.,  Koomen,  M.J.,  Korendyke,  C.M.,  Michels,  D.J.,  Moses, 
J.D.,  Socker,  D.G.,  Dere,  K.P.,  Lamy,  P.L.,  Llebaria,  A.,  Bout,  M.V.,  Schwenn,  R.,  Sim- 
nett,  G.M.,  Bedford,  D.K.,  Eyles,  C.J.:  1995,  The  large  angle  spectroscopic  coronagraph 
(LASCO).  Solar  Phys.  162(1-2),  357-402.  doi:10.1007/BF00733434. 

Cremades,  H.,  Bothmer,  V.,  Tripathi,  D.:  2006,  Properties  of  structured  coronal  mass  ejections 
in  solar  cycle  23.  Adv.  Space  Res.  38,  461-465.  doi:10.1016/j.asr.2005.01.095. 

Gosling,  J.T.:  1990,  Geophysical  monograph.  In:  Priest,  E.R.,  Lee,  L.C.,  Russel,  C.T.  (eds.) 
Coronal  mass  ejections  and  magnetic  flux  ropes  in  interplanetary  space  58,  343-364. 

Howard,  R.A.,  Moses,  J.D.,  Vourlidas,  A.,  Newmark,  J.S.,  Socker,  D.G.,  Plunkett,  S.P.,  Ko¬ 
rendyke,  C.M.,  Cook,  J.W.,  Hurley,  A.,  Davila,  J.M.,  Thompson,  W.T.,  St  Cyr,  O.C., 
Mentzell,  E.,  Mehalick,  K.,  Lemen,  J.R.,  Wuelser,  J.P.,  Duncan,  D.W.,  Tarbell,  T.D., 
Wolfson,  C.J.,  Moore,  A.,  Harrison,  R.A.,  Waltham,  N.R.,  Lang,  J.,  Davis,  C.J.,  Eyles, 
C.J.,  Mapson-Menard,  H.,  Simnett,  G.M.,  Halain,  J.P.,  Defise,  J.M.,  Mazy,  E.,  Rochus,  P., 
Mercier,  R.,  Ravet,  M.F.,  Delmotte,  F.,  Auchere,  F.,  Delaboudiniere,  J.P.,  Bothmer,  V., 
Deutsch,  W.,  Wang,  D.,  Rich,  N.,  Cooper,  S.,  Stephens,  V.,  Maahs,  G.,  Baugh,  R.,  Mc- 
Mullin,  D.,  Carter,  T.:  2008,  Sun  Earth  connection  coronal  and  heliospheric  investigation 
(SECCHI).  Space  Sci.  Rev.  136,  67-115.  doi:10.1007/sll214-008-9341-4. 

Hu,  Q.,  Sonnerup,  B.U.O.:  2002,  Reconstruction  of  magnetic  clouds  in  the  solar  wind: 
Orientations  and  configurations.  J.  Geophys.  Res.  107,  1142.  doi:10.1029/2001JA000293. 

Huttunen,  K.E.J.,  Koskinen,  H.E.J.,  Schwenn,  R.:  2002,  Variability  of  magnetospheric  storms 
driven  by  different  solar  wind  perturbations.  J.  Geophys.  Res.  107(A7),  1121-1128. 
doi:10.1029/2001JA900171 . 

Isavnin,  A.,  Kilpua,  E.K.J.,  Koskinen,  H.E.J.:  2011,  Grad-Shafranov  reconstruction 
of  magnetic  clouds:  overview  and  improvements.  Solar  Phys.  273(1),  205-219. 
doi:  10. 1007 /si  1207-01 1-9845-z . 

Jian,  L.,  Russel,  C.T.,  Luhmann,  J.G.,  Skoug,  R.M.:  2006,  Properties  of  interplane¬ 
tary  coronal  mass  ejections  at  one  AU  during  1995-2004.  Solar  Phys.  239,  393-436. 
doi:  10. 1007 /sll207-006-0133-2 . 

Kaiser,  M.L.,  Kucera,  T.A.,  Davila,  J.M.,  Cyr,  O.C.S.,  Guhathakurta,  M.,  Crhistian, 
E.:  2008,  The  STEREO  mission:  An  introduction.  Space  Sci.  Rev.  136,  5-16. 
doi:10.1007/sll214-007-9277-0. 

Kilpua,  E.K.J.,  Pomoell,  J.,  Vourlidas,  A.,  Vainio,  R.,  Luhmann,  J.,  Li,  Y.,  Schroeder,  P., 
Galvin,  A.B.,  Simunac,  K.:  2009,  STEREO  observations  of  interplanetary  coronal  mass 
ejections  and  prominence  deflection  during  solar  minimum  period.  Annales  Geophysicae 
27,  4491-4503.  doi:10.5194/angeo-27-4491-2009. 

Krall,  J.:  2007,  Are  all  coronal  mass  ejections  hollow  flux  ropes?  Astrophys.  J.  657(1),  559- 
566.  doi:10. 1086/510191. 

Lynch,  B.J.,  Antiochos,  S.K.,  Li,  Y.,  Luhmann,  J.G.,  DeVore,  C.R.:  2009,  Rota¬ 
tion  of  coronal  mass  ejections  during  eruption.  Astrophys.  J.  697,  1918-1927. 
doi:  10. 1088/0004-637X/697/2/1918 . 

MacQueen,  R.M.,  Hundhausen,  A.J.,  Conover,  C.W.:  1986,  The  propagation  of  coronal  mass 
ejection  transients.  J.  Geophys.  Res.  91(A1),  31-38.  doi:10.1029/JA091iA01p00031. 

Nieves-Chinchilla,  T.,  Colaninno,  R.,  Vourlidas,  A.,  Szabo,  A.:  2012,.  J.  Geophys.  Res.  in 
press. 

Plunkett,  S.P.,  Thompson,  B.J.,  Cyr,  O.C.S.,  Howard,  R.A.:  2001,  Solar  source  regions  of 
coronal  mass  ejections  and  their  geomagnetic  effects.  J.  Atmos.  Solar  Terr.  Phys.  63, 
389-402.  doi:10.1016/S1364-6826(00)00166-8. 


SOLA:  Isavnin_etal_2012a.tex;  5  December  2012;  4:31;  p.  13 


A.  Isavnin,  et  al. 


Rodriguez,  L.,  Mierla,  M.,  Zhukov,  A.N.,  West,  M.,  Kilpua,  E.:  2011,  Linking  of  remote- 
sensing  and  in  situ  observations  of  CMEs  using  STEREO.  Solar  Phys.  270,  561-573. 
doi:  10. 1007/sll207-01 1-9784-8. 

Thernisien,  A.:  2011,  Implementation  of  the  graduated  cylindrical  shell  model  for  the  three- 
dimensional  reconstruction  of  coronal  mass  ejections.  Astrophys.  J.  Suppl.  Ser.  194, 
33(6pp).  doi:  10. 1088 / 0067-0049 / 194/2/33 . 

Thernisien,  A.,  Vourlidas,  A.,  Howard,  R.A.:  2009,  Forward  modeling  of  coronal  mass  ejections 
using  STEREO/SECCHI  data.  Solar  Phys.  256,  111-130.  doi:10.1007/sll207-009-9346-5. 

Thompson,  W.T.:  2006,  Coordinate  systems  for  solar  image  data.  Astronomy  and  Astroph. 
449(2),  791-803.  doi:10. 1051/0004-6361:20054262. 

Torok,  T.,  Kliem,  B.:  2003,  The  evolution  of  twisting  coronal  magnetic  flux  tubes,  a  406, 
1043-1059.  doi:10. 1051/0004-6361:20030692. 

Tsurutani,  B.T.,  Gonzalez,  W.D.,  Tang,  F.,  Akasofu,  S.I.,  Smith,  E.J.:  1988,  Origin  of  inter¬ 
planetary  southward  magnetic  fields  responsible  for  major  magnetic  storms  near  solar  max¬ 
imum  (1978-1979).  J.  Geophys.  Res.  93(A8),  8519-8531.  doi:10.1029/JA093iA08p08519. 

Vourlidas,  A.,  Colaninno,  R.,  Noeves-Chinchilla,  T.,  Stenborg,  G.:  2011,  The  first  observation 
of  a  rapidly  rotating  coronal  mass  ejection  in  the  middle  corona.  Astrophys.  J.  Lett.  733, 
L23.  doi:10.1088/2041-8205/733/2/L23. 

Wang,  Y.M.,  Ye,  P.Z.,  Wang,  S.,  Zhou,  G.P.,  Wang,  J.X.:  2002,  A  statistical  study  on  the 
geoeffectiveness  of  Earth-directed  coronal  mass  ejections  from  March  1997  to  December 
2000.  J.  Geophys.  Res.  107(A11),  1340.  doi:  10. 1029/2002 JA009244. 

Wang,  Y.,  Shen,  C.,  Wang,  S.,  Ye,  P.:  2004,  Deflection  of  coronal  mass  ejection  in  the  inter¬ 
planetary  medium.  Solar  Phys.  222,  329-343.  doi:10. 1023/B:SOLA. 0000043576. 21942. aa . 

Wang,  Y.,  Chen,  C.,  Gui,  B.,  Shen,  C.,  Ye,  P.,  Wang,  S.:  2011,  Statistical  study  of  coronal 
mass  ejection  source  locations:  Understanding  CMEs  viewed  in  coronagraphs.  J.  Geophys. 
Res.  116(A04104).  doi:10.1029/2010JA016101. 

Yurchyshyn,  V.:  2008,  Relationship  between  EIT  post  eruption  arcades,  coronal  mass 
ejections,  coronal  neutral  line  and  magnetic  clouds.  Astrophys.  J.  Lett.  675,  L49. 
doi:  10. 1086/533413. 

Yurchyshyn,  V.,  Abramenko,  V.,  Tripathi,  D.:  2009,  Rotation  of  white-light  coronal  mass 
ejection  structures  as  inferred  from  LASCO  coronagraph.  Astrophys.  J.  705,  426-435. 
doi:  10. 1088 /0004-637X/705/1/426 . 

Zhang,  J.,  Richardson,  I.G.,  Webb,  D.F.,  Gopalswamy,  N.,  Huttunen,  E.,  Kasper,  J.C.,  Nitta, 
N.V.,  Poomvises,  W.,  Thompson,  B.J.,  Wu,  C.-C.,  Yashiro,  S.,  Zhukov,  A.N.:  2007,  Solar 
and  interplanetary  sources  of  major  geomagnetic  storms  (Dst  100  nT)  during  19962005. 
J.  Geophys.  Res.  112(A10102).  doi:10.1029/2007JA012321. 

Zurbuchen,  T.H.,  Richardson,  I.G.:  2006,  In-situ  solar  wind  and  magnetic  field  sig¬ 
natures  of  interplanetary  coronal  mass  ejections.  Space  Sci.  Rev.  123,  31-43. 
doi:10.1007/sll214-006-9010-4. 


SOLA:  Isavnin_etal_2012a.tex;  5  December  2012;  4:31;  p.  14 


